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ABSTRACT 

Aims. A variation of the fundamental constants is expected to affect the thermonuclear rates important for stellar nucleosynthesis. In 
particular, because of the very small resonant energies of 8 Be and l2 C, the triple a process is extremely sensitive to any such variations. 
Methods. Using a microscopic model for these nuclei, we derive the sensitivity of the Hoyle state to the nucleon-nucleon potential 
allowing for a change in the magnitude of the nuclear interaction. We follow the evolution of 15 and 60 Mq zero metallicity stellar 
models, up to the end of core helium burning. These stars are assumed to be representative of the first, Population III stars. 
Results. We derive limits on the variation of the magnitude of the nuclear interaction and model dependent limits on the variation of 
the fine structure constant based on the calculated oxygen and carbon abundances resulting from helium burning. The requirement 
that some 12 C and ls O be present are the end of the helium burning phase allows for permille limits on the change of the nuclear 
interaction and limits of order 10~ 5 on the fine structure constant relevant at a cosmological redshift of z ~ 15 - 20. 
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1. Introduction 

The equivalence principle is a cornerstone of metric the ories o f 
gravitation and in particular of general relativity d"Willl 1 1 993b . 
This principle, including the universality of free fall, the local 
position and Lorentz invariances, postulates that the local laws 
of physics, and in particular the values of the dimensionless con- 
stants such as the fine structure constant a em = e 2 /4nsofic, must 
remain fixed, and thus be the same at any time and in any place. 
It follows that by testing the constancy of fundamental constants 
one actually performs a test of General Relativity, that can be 
ext ended on astrophy sical and cosmological scales (for a review, 
see lUzanll200ll2009ah 

We define a fundamental const ant as any free parameter of 
the fundamental theories at hand d Weinberg [ljii3t iDuffl 120021: 
iDuff et all 12001 iBarrowl 120021: lUzan & Leclercqll2008l) . These 
parameters are contingent quantities that can only be measured 
and are assumed constant since (i) in the theoretical framework 
in which they appear, there is no equation of motion for them and 
they cannot be deduced from other constants and (ii) if the the- 
ories in which they appear have been validated experimentally, 
it means that, at the precision of the experiments, these parame- 
ters have indeed been checked to be constant. Hence, by testing 
for their constancy we extend our knowledge of the domain of 
validity of the theories in which they appear. In that respect, as- 
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trophysics and cosmology allow one to probe larger time-scales, 
typically of the order of the age of the universe. 

One can, however, question the constancy of these dimen- 
sionless numbers and the physics which determines their value. 
This s ends us back to the phenomenological argument by Dirac 
d!937l) . known as the 'Large Number Hypothesis', according to 
which the dimensionless ratio Gm e m p /hc, or simply G in atomic 
units, should decreas e as th e inverse of the age of the universe, 
followed by iJordanl ( 1 19371) who formulated a field theory in 
which both the fine structure constant and the gravitational con- 
sta nt wer e repla ced by dynamical fields. It was soon pointed out 
by iFierj dl956l) that astronomical observations can set strong 
constraints on the variations of these constants. This paved the 
way to two complementary directions of research on the funda- 
mental constants. 

On the one hand, from a theoretical perspective, many 
theories involving "varying constants" have been designed. 
This is in particular the case of theories involvin g extra- 
dimensions, such as the Kaluza-Klein mechanism dKaluzal 
11921b iKleinl [T926) and string theory, in which all the con- 
stants (including gauge, Yukawa and gravitational couplings) 
are dynamical quantities ( )Wu & Wand 1 19861 IWetterichl ll988: 
Taylo r & Venezianoll 19881: [Witten 1984), or in theories such as 
scalar- tensor theories of gravity ( Jordanl 119491 : IBrans & Dickd 
1 19611: iDamour & E sposito-Farese 199 2[) and in m a ny models 
of qu i ntessence dUzanlll999l: IDamour et al, 2002a b; [Wetterichl 
2003; iLee et alJ 120041: iRiazuelo & Uzanl I2001 T aiming at ex- 
plaining the acceleration of the universe by the dynamics 



2 



Ekstrom et al.: Effect of the variation of fundamental constants on Pop III star evolution 



of a scalar field. It is impingent on these models to ex- 
plain why the constants are so c onstant today and provide a 
mechanism for fixing their value (Damour & Nordtvedtl ll993; 
Dam our & Polyakovll 19941) . In this respect, testing for the con- 
stancy of the fundamental constants is one of the few windows 
on these theories. 

On the other hand, from an experimental and observational 
perspective, the variation of various constants have been 
severely constrained. This is the case of the fine structure con- 
stant for which the constraint a em /a em = (-1.6+2.3)xl0~ 17 yr~ 1 
at z — has been obtained from the c omparison of aluminium 
and mercury single-ion optical clocks (Rosenba nd et al.l l2008). 
Over a longer timescale, it was demonstrated that a em cannot 
have varied by mo re than 10~ 7 over the last 2 Gyr from the 
Oklo phenomenon dShlvakhterlll976t iDamour & Dvsonlll996t 
iFuiii et alj|2000t lOlive et al J 120021: iPetrov et alj|200a and over 
the la s t 4.5 Gyr from meteor i te dating (Dickd [19591: iDvsonl 
119721: IFuiii & Iwamotol l2003t lOlive et alj 120041) . At higher 
redshift, 0.4 < z < 3.5, there are conflicting reports of a an 
observed variation of a em from qua sar absorpt i on sys tems. 
Using the m a ny-mu l tiplet method, Webb et all (120011) and 
Mur phy et al.l (120031 120071) claim a statistically significant 
variation Aa em /a em - (-0.54 + 0.12) x 10~ 5 , indicating a 
smaller value of a em in the past. More recent observations 
taken at VLT/UVES using the many multiplet method have 
not been able to duplicate the previous resu l t dChand et al. | 
2004; ISrianand et ail 12004 lOuast et al.l 120041: ISrianand et al.l 
2007). The use of Fe lines in lOuast et all d2004l) on a sin- 
gle absorber found Aa em /a em = (-0.05 ± 0.17) x 10" 5 . 
However, since the previous result relied on a statistical 
average of over 100 absorbers, it is not clea r that these two 
results are in contradiction. In IChand et al.l d2004l) . the use 
of Mg and Fe lines in a set of 23 systems yielded the result 
Aa em /a em = (0.01 +0.15)x 10~ 5 and therefore represents a more 
significant disagreement and can be used to set very stringent 
limits on the possible variation in a em . A pu rely astrophysical 
explan ation for these results is also possible dAshenfelter et alj 
2004ab). At lar ger redshifts, constraints at the percent level 
have been obtained from the observation of the temperature 
anisot ropies of cosmic m i crowave background at (z ~ 10 3 ) 
iMartins et all 120041: Istefanescul l2007t iNakashima et alj 
IScoccola et al.ll2008l) and from big bang nucleosynthesis 
(BBN ) ( z ~ 10 10 ) (a 
19951: iBergstrom et; 
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Muller et all 12004 Ichikawa & Kawasaki 2004; lLandau et alJ 
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2006; Coc et al. 2007; D ent et al.l2007l) . We refer to ltfzanl d2003l 



2009b) for recent reviews on this topic. For the time being, 
there is no constraint on a em for redshifts ranging from 4 to 
10 3 although it has been proposed that 21 cm observations may 
allow one to fill in the range 30 < z < 100 dKhatri & Wa ndelt 
l2007h . 

This article focuses on the effect of the possible variation 
of the fundamental constants on the stellar evolution of early 
stars, hence possibly providing constraints in a domain of red- 
shifts where no such c onstraint is avai lable. A similar issue was 
ac tually consider ed bv lGamowl d 19671) (see also the recent work 
bv lAdamsl d2008l) who showed that the evolution of the Sun was 
able to exclude the Dirac model of a varying gravitational con- 
stant. In this case, non-gravitational physics is kept unchanged 
and the evolution of the star is affected only by the modifica- 
tion of gravity. Changing the non-gravitational sector has more 
drastic implications on stellar physics since the nuclear physics 



and thus the cross-sections and reaction rates of all the processes 
sh ould be modified. 

iRozentail dl988l) argued that the synthesis of complex ele- 
ments in stars (mainly the possibility of the Sa-reaction as the 
origin of the production of 12 C) sets constraints on the values 
of the fine structure and strong coupling constants. There have 
been several studies on the s ensitivity of carbon production to 
the underlying nu clear rates ( Livio et al.lll989t lFairbairn|[ l999: 
Csoto et all 1200 ll lO berhummer et al. 2000{ lOberhummer et all 
l2003b ISchlattl et alj 12004 iTur et al] 120071) . The production of 
12 C in stars requires a triple tuning: (i) the decay lifetime of 
8 Be, of order 10~ 16 s, is four orders of magnitude longer than 
the time for two a particles to scatter, (ii) an excited state of the 
carbon lies just above the energy of 8 Be + a and finally (iii) the 
energy level of 16 at 7. 1 197 MeV is non resonant and below the 
energy of 12 C + a, at 7.1616 MeV, which ensures that most of 
the carbon synthesized is not destroyed by the capture of an a- 
particle. The existence of t his excited state of 12 C was actually 
predicted by Hovld d!95 4) and then observed at th e predicted 
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119571) . The variation of any constant which would modify the 
energy of this resonance, known as the Hoyle level, would dra- 
matically affect the production of carbon. 

Qualitatively, and perhaps counter-intuitively, if the energy 
level of the Hoyle level were increased, 12 C would probably be 
rapidly processed to 16 since the star would, in fact, need to 
be hotter for the triple-a reaction to be triggered. On the other 
hand, if it is decreased very little oxygen will be produced. From 
the general expression of the reaction rate (see below for details, 
definitions of all the quantities entering this expression, and a 
more accurate computation) 



A 3a = 3 3/ X 3 



M a k B rj h exp 



Qa 



where Q aaa ~ 380 keV is the energy of the resonance, one 
deduces that the sensitivity of the reaction rate to a variation of 

Qaaa IS 

dhU 3 „ Q aaa (-4.4] 



din Q a 



T 9 



where Tg = r/ 10 9 K. This effect w a s inves tig ated by Csoto et al. 
(200lh and lOberhummer et al] (120001) ; lOberhummer et al. 
(120031) who related the variation of Q aaa to a variation of the 
strength of the nucleon-nucleon (N-N) interaction. Focusing on 
the C/O ratio in red giant stars (1.3, 5 and 20 Mq with so- 
lar metallicity) up t o thermally pulsing asymptotic giant branch 
stars (TP-AGB) (lOberhummer et al] l2000t lOberhummer et al.l 
120031) and in low, intermediate and high mass stars (1. 3, 5, 15 
and 25 M Q with solar metallicity) dSchlattl et al]|2004 . it was 
estimated that outside a window of 0.5% and 4% for the values 
of the strong and electromagnetic forces respectively, the stellar 
production of c arbon or oxygen wi ll be reduced by a factor 30 to 
1000 (see also lPochet etail d!991l) ). 

Indeed, modifying the energy of the resonance alone is not 
realistic since all cross-sections, reaction rates and binding en- 
ergies etc. should be affected by the variation of the constants. 
One could indeed have started by assuming independent varia- 
tions of all these quantities but it is more realistic (and hence 
more model-dependent) to try to deduce their variation from a 
microscopic model. Our analysis can then be outlined in three 
main steps: 

1. Relating the nuclear parameters to fundamental constants 
such as the Yukawa and gauge couplings, and the Higgs 
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vacuum expectation value. This is a difficult step because 
of the intricate structure of QCD and its role in low energy 
nuclear reactions, as in the case of BBN. The nuclear param- 
eters include the set of relevant energy levels (including the 
ground states), binding energies of each nucleus and the par- 
tial width of each nuclear reaction. This involves a nuclear 
physics model of the relevant nuclei (mainly 4 He, 8 Be, 12 C, 
and 16 for our study). 

2. Relating the reaction rates to the nuclear parameters, which 
implies an integration over energy of the cross-sections. 

3. Deducing the change in the stellar evolution (lifetime of the 
star, abundance of the nuclei, Hertzprung-Russel (H-R) dia- 
gram, etc.). This involves a stellar model. 

Let us summarize the main hypothesis of our work for each 
of these steps. 

The first step is probably the most difficult. We shall adopt a 
phenomenological description of the different nuclei based on a 
cluster model in which the wave functions of the 8 Be and 12 C nu- 
clei are approximated by a cluster of respectively two and three a 
wave functions. When solving the associated Schrodinger equa- 
tion, we will modify the strength of the electromagnetic and nu- 
clear N-N interaction potentials respectively by a factor ( 1 + S a ) 
and (1 + 6nn) where 5 a and 6nn are two small dimensionless 
parameters that encode the variation of the fine structure con- 
stant and other fundamental couplings. At this stage, the rela- 
tion between 6nn an d the gauge and Yukawa couplings is not 
known. This will allow us to obtain the energy levels, includ- 
ing the binding energy, of 2 H, 4 He, 8 Be, 12 C and the first J" = 
+ 12 C excited energy level. Note that all of the relevant nuclear 
states are assumed to be interacting alpha clusters. In a first ap- 
proximation, the variation of the a particle mass cancels out. The 
partial widths (and lifetimes) of these states are scaled from their 
experimental laboratory values, according to their energy depen- 
dence. 6nn is used as a free parameter. The dependence of the 
deuterium binding energy on 6nn then offers us the possibility of 
relating this parameter to the gauge and Yukawa couplings if one 
matches this prediction to a potential model via the cr and u> me- 
son masse s (Flamb aum & Shurvakl 120031; iDmitriev et al.| [2004; 
ICoc et al l 120071: iDamour & DonoghueT 2008 ) or the pion mass , 
as suggested bylYoo & Scherrerl (12003): Epelbau m et al.l(l2003l) : 
lBean& Savage! (120031) . 

The second step requires an integration over energy to de- 
duce the reaction rates as functions of the temperature and of the 
new parameters 8 a and 6mn- 

The third step involves stellar models and in particular some 
choices about the masses and initial metallicity of the stars. 
While theoretically uncertain, it is usually thought that the first 
stars were massive; how ever, their mass ran ge is presently un- 
known, (for a review, see lBromm e t al. (2009)). In a hierarchical 
scenario of structure formation, they were formed a few x 10 8 
years after the big bang, that is at a redshift of z ~ 10-15 
with zero metallicity (so that we can use the BBN abundances as 
initial conditions, see Sect. @}. We thus focus on Population III 
stars (Pop III) with typical masses, 15 and 60 M@, assuming no 
rotation. Our computation is stopped at the end of core helium 
burning. 

The final step would be to use these predictions to set con- 
straints on the fundamental constants, using stellar constraints 
such the C/O ratio which is in fact observable in very metal poor 
stars. While this article can be seen as a theoretical investigation 
that describes the expected effect of a variation of the funda- 
mental constants, it also sheds some interesting light on stellar 
physics and its sensitivity to fundamental physics. 



The article is logically organized as follows. Section 2 re- 
calls the basis of the Sa-reaction, Section 3 describes the nuclear 
physics modeling (first step), Section 4 is devoted to stellar im- 
plications and Section 5 to the discussion. Technical details are 
gathered in the appendices. 

2. Stellar carbon production 

The 3a process is one of the most delicate of all reactions in nu- 
clear astrophysics. It is also one of the most influential since it 
bypasses the deep gap between BBN and stellar nucleosynthesis. 
More specifically, BBN stops at mass 7 ( 7 Li) due to the lethal 
instability of parent nuclei with strongly bound 4 He offsprings, 
namely nuclei with masses 5 and 8. The 3a reaction allows nu- 
clear complexity to proceed up to uranium through core collapse 
supernova explosions. Note that the gap between BBN and stel- 
lar nucleosynthesis is filled by non-thermal processes (spallative 
processes induced by galactic cosmic rays) producing 6,7 Li, 9 Be 
and 101I B. Once these nuclear obstacles are overcome, the phys- 
ical conditions within stars allow nuclear production of elements 
from carbon and beyond. As such, the 3a reaction is the first step 
of helium burning which is followed in massive stars by C, Ne, 
O, and Si burning and then explosive nucleosynthesis. The mass 
of the C-O core and the C/O ratio at the end of helium burning 
is important for determining i) the subsequent phases of stellar 
evolution and nucleosynthesis since it fixes the mass of the iron 
core and ii) the final fate of stars (black holes, neutron stars, or 
white dwarfs). In particular for the first stars, the 3a reaction is 
of great importance since no metals have yet been formed and 
the CNO cycle cannot proceed as usual. Unfortunately, the 3a 
reaction is a two step sequential process and the & Be(a,y) l2 C 
cross section has not been measured directly in the laboratory. 
Indeed, the 8 Be lifetime (about 10~ 16 sec) is so short that such a 
measurement is not currently feasible. 

Consequently, the C and O abundances at the end of helium 
burning is very sensitive to small variations of the 3a reaction 
rate. In this context, any anomalous abundance of C and O in 
very metal poor stars could potentially be taken as an indication 
of the variation of the nucleon - nucleon interaction and there- 
fore of either or both of the electromagnetic and strong coupling 
constants. 

In our analysis, we focus on the C/O ratio. It is of interest, 
therefore, to comment on the destruction of carbon (production 
of oxygen) as well as the destruction of oxygen. If the reaction 
following the 3a process, namely 12 C(or, y) I6 0, is sufficiently 
fast, then most a particles would be converted to ie O or heavier 
nuclei with little I2 C left at the end of helium burning. However, 
the fact that in general the C/O ratio in the Universe is about 
0.4 suggests that the 12 C(a, y) lfi O reaction is sufficiently slow 
that some 12 C remains after helium exhaustion. The presence of 
comparable quantities of C and O implies also that the subse- 
quent lfi O( a, y) 20 Ne reaction is not too fast , otherwise O would 
be converted to Ne or heavier nuclei and little O would survive 
during helium burning. We would like to stress the importance 
of the nuclear balance between C and O. The observation of C/O 
in very metal poor stars may hold the key to any variation in the 
chain of processes described above. 

In Fig. Q] we show the low energy level schemes of the nu- 
clei participating to the 4 He(aa, y) I2 C reaction: 4 He, 8 Be and 
I2 C. The triple a process begins when two alpha particles fuse 
to produce a 8 Be nucleus whose lifetime is only ~ 10~ 16 s but 
is sufficiently long so as to allow a second alpha capture into 
the second excited level of 12 C, at 7.65 MeV above the ground 
state (of 12 C). In the following, we shall refer to the successive 
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Fig. 1. Level scheme showing the key levels in the triple a pro- 
cess. 



a captures as first and second steps, that is aa <-> Be+y and 
8 Be+a <-> 12 C* — » 12 C+y. The excited state of 12 C c orresponds 
to an I — resonance, as postulated bv lHovldd 19541) in order to 
increase the cross section during the helium burning phase. This 
level decays to the first excited level of 12 C at 4.44 MeV through 
an E2 (i.e. electric with 1-2 multipolarity) radiative transition 
as the transition to the ground state (0| — > OJ) is suppressed 
(pair emission only). At temperatures above Tg =s 2, which are 
not relevant for our analysis and therefore not treated, one should 
also consider other possible levels above the a threshold. 
We define the following energies: 

- £«( 8 Be) as the energy of the 8 Be ground state with respect 
to the a+a threshold; 

- Er{ C) as the energy of the Hoyle level with respect to the 
8 Be+a threshold, i.e. E R ( 12 C) = 12 C(0 2 + )+2 t? ( 12 C) where 
12 C(02 + ) is the excitation energy and £2u-( 12 C) is the a par- 
ticle separation energy; 

- Qaaa as the energy of the Hoyle level with respect to the 3a 
threshold so that 



Q aaa = £ fi ( 8 Be) + £ S ( 12 C). 



(1) 



- r a ( 8 Be) as the partial width of the beryllium decay 
{aa <-» 8 Be+y); 

- r rit? ( 12 C) as the partial widths of & Be+a ^ 12 C* -+ 12 C+y. 

Their standard values are given in Table [1] 

Assuming i) thermal equilibrium between the 4 He and 8 Be 
nuclei, so that their abundances are related by the Saha equation 
and ii) the sharp resonance approximation for the alpha capture 
on 8 Be, the 4 He(aa , y) 12 C rate can be expressed (Nomoto et al. 
1 1 985b Illiadisll2007h as: 



N 2 A {a-v) aaa 



3 i/2 6Nl 



In 



M a k B T 



*s i-Q 

n wyexp 



k B T 



(2) 



Table 1. Nuclear data for the two steps o f the 3a-reaction . 
See text for the definitions of the quantit ies dTillev et alj [2004; 
Aizenberg-Selove 1990: [Audi et allhoolh . 



nucleus 


J" 


E R (keV) 


T„ (eV) 


T r (meV) 


"Be 


+ 


91.84±0.04 


5.57+0.25 




12 C 


a 


287.6+0.2 


8.3 ± 1.0 


3.7 ±0.5 



withw = 1 (spin factor), y = r r ( 12 C)r ff ( 12 C)/(r r ( I2 C)+r a ,( 12 C)) 
« T r ( 12 C) for present day values, and M a is the mass of the a 
nucleus. 

During he lium burning , the only other important reaction is 
12 C(a,y) 16 dlliadisl 120071) which transforms 12 C into lfl O. Its 
competition with the 3a reaction governs the I2 C/ ie O abundance 
ratio at the end of the helium burning phase. Even though, the 
precise value of the 12 C(or, y) 16 S-factoiQ is still a matter of 
debate as it relies on an extrapolation of experimental data down 
to the astrophysical energy (^300 keV), its energy dependence 
is much weaker than fthat of the 3a reaction. Indeed, as it is 
dominated by broad resonances, a shift of a few hundred keV in 
energy results in a S-factor variation of much less than an order 
of magnitude. For this reason, we can safely neglect the effect 
of the l2 C(a, y) 16 reaction rate variation when compared to the 
variation in the 3a rate. Similar considerations apply to the rate 
for I6 O(a,y) 20 Ne. 

During hydrogen burning, the pace of the CNO cycle is given 
by the slowest reaction, 14 N(p,y) 15 0. Its 5-factor exhibits a well 
known resonance at 260 keV which is normally outside of the 
Gamow energy window (s;100 keV) but a variation of the N-N 
potential could shift its position downward, resulting in a higher 
reaction rate and more efficient CNO H- burning. 

3. Microscopic determination of the 3a rate 

3. 1 . Description of the cluster model 

In order to analyze the sensitivity of the triple a reac- 
tion to a variation of the strength of the electromagnetic 
and NN interactions, w e use a microscopic model (see 
IWildermufh & Tangl 11977b iKorennov & Descouvemontl 120041 
and references therein). In such an approach, the wave function 
of a nucleus with atomic number A, spin J, and total parity n is a 
solution of a Schrodinger equation with a Hamiltonian given by 



h =Tj Tird + Z y(r " ) - 

i !>J=1 



(3) 



r(r,) is the kinetic energy of nucleon i. The nucleon-nucleon 
interaction V(r (/ ) depends only on the set of relative distances 



V(ry) 



- . j. It can be decomposed as 
Vc(ry) + V N (r u ), 



(4) 



where the potential Vc( r ) arises from the electromagnetic inter- 
action and Vjv(r) from the nuclear interaction. The expression 
for Vn is detailed in Appendix[A] The eigenstates x ^ JMn with en- 
ergy E ln of the system are solutions, as usual, of the Schrodinger 
equation associated with the Hamiltonian given in Eq. (0, 



JMk 



Jn\r<JMK 



(5) 



1 The astrophysical 5-factor is just the cross section corrected for 
the effect of the penetrability of the Coulomb barrier and other trivial 
effects. 
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The total wave function m JMn is a function of the A - I coordi- 
nates Yh, 

When A > 4, no exact solutions of Eq. (O can be found and 
approximate solutions have to be constructed. For those cases, 
we use a cluster approximation in which x i> JMn is written in 
terms of a-nucleus wave functions. Because the binding energy 
of the a particle is large, this approach has been shown to be 
well adapted to c l uster states, and in particular to 8 Be and 12 C 
(Kamimural ll98ll : ISuzuki et al.ll2.Q08l) . In the particular case of 
these two nuclei, the wave functions are respectively expressed 
as 

v{£ = Mt> a <t> a g{ M "(p) 

Vif = Jlcp^^g^ifi, R), (6) 

where <f> a is the a wave function, defined in the 0s shell model 
with an oscillator parameter b; J[ is the antisymmetrization op- 
erator between the A nucleons of the system. For two-cluster 
systems, the wave function gi (p) depends on the relative coor- 
dinate p between the two a particles. For three-cluster systems, 
R is the relative distance between two a particles, and p the rel- 
ative coordinate between the third a particle and the 8 Be center 
of mass. The relative wave functions, g2 and g$, are obtained by 
solving the Schrodinger equation (0. 

One then needs to specify the nucleon-nucleon poten- 
tial yjy(r i; ). We shall use the microscopic interaction model 
dThompson et al.l 19771) which contains one linear parameter (ad- 
mixture parameter u), whose standard value is u = 1. It can 
be slightly modified to reproduce important inputs, such as the 
resonance energy of the Hoyle state. The binding energies of 
the deuteron (-2.22 MeV) and of the a particle (-24.28 MeV) 
do not depend on u. For the deuteron, the Schrodinger equa- 
tion is solved exactly. More details about the model are given in 
AppendixlAl 

To take into account the variation of the fundamental con- 
stants, we introduce the parameters 5 a and 5nn to characterize 
the change of the strength of the electromagnetic and nucleon- 
nucleon interaction respectively. This is implemented by modi- 
fying the interaction potential (|4]i so that 

V(r y ) = ( 1 + S a ) V c (r y ) + ( 1 + 6 NN )V N (r u ). (7) 

Such a modification will affect Bjy and the energy levels of 
8 Be and 12 C simultaneously. Of course, one could have imag- 
ined a more complex parameterization of the variations (e.g., 
by varying all quantities in Eq. ( 1A2I )). but since we expect to 
consider only small variations of all quantities, as an approx- 
imation, the system should be linear in the variations and our 
approach should be sufficient for extracting the physical effects 
of any such small variation. 

3.2. Sensitivity of the nuclear parameters 

For each set of values (5 a , <5nn) we solve Eq. (0 with the inter- 
action potential (Q. We emphasize that the parameter u is de- 
termined from the experimental 8 Be and 12 C(0£) energies (u = 
0.954). We assume that 6 NN varies in the range [-0.015, 0.015]. 

First, concerning the deuteron, this analysis implies that its 
binding energy scales as 

ABd/Bd = 5.716 x SffN . (8) 

Second, concerning 8 Be and 12 C, we can extract the sensitivity 
of E R ( & Be) and E R ( l2 C). They scale as 

£ R ( 8 Be) = (0.09184 - 12.208 x 6 NN ) MeV (9) 



ABj/B D 

-0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 




-0.015 -0.01 -0.005 0.005 0.01 0.015 



Fig. 2. Variation of the resonance energies as a function of 6nn- 
The symbols represent the results of the microscopic calculation 
while the lines correspond to the adopted linear relationship be- 
tween Er and Snn- 

and 

E R ( 12 C) = (0.2876 - 20.412 x S NN ) MeV. (10) 

The numerical results for the sensitivities of £s( 8 Be) and 
£r( 12 C) to 6nn as well as the above linear fits are shown in 
Fig. [2] The effect of 6 a on these quantities is negligible. Note 
that for 6nn Z 0.007, £ , «( 8 Be) is negative and 8 Be would become 
stable. Using the bijective relation ([8]l between Bo and Snn we 
can also express our results as 

£«( 8 Be) = (0.09184 - 2.136 x ABo/Bo) MeV, (11) 
£ S ( 12 Q = (0.2876 - 3.570 x ABo/Bp) MeV. (12) 

It follows that the energy of the Hoyle level with respect to the 
3a threshold (and not with respect to 8 Be+a threshold) is given 
by (see Eq. ([TJ) 

Qaaa = (0.37945 - 5.706 x ABd/Bd) MeV, (13) 
= (0.37945 - 32.620 x 6 NN ) MeV. (14) 

To estimate the effect of 6 a in Eq. (7), we can approx- 
imate the Coulomb energy by (3/5)Z(Z - l)a em hc/R c where 
R c = 1.3A 1/3 fm which gives 9 MeV for 12 C and 0.9 MeV 
for 4 He. The variation of Q aaa is thus of the order of +6 MeV 
x 5 a . The direct effect of S a is thus of opposite sign but con- 
siderably less important. This is in qualitative agreeme nt with 
lOberhummer et all (120001) and lOberhummer et al.T(l2001l) . 

3.3. Sensitivity of the 3a-reaction rate 

The method described above provides a consistent way to eval- 
uate the sensitivity of the 3a-reaction rate to a variation of the 
co nstants. This rate has been computed numerically as explained 
in lAngulo et al.l d 19991) and as described in Appendix [B] where 
both an analytical approximation valid for sharp resonances and 
a numerical integration are performed. 
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Fig. 3. The ratio between the 3a rate obtained for -0.009 < 
Snn ^ +0.006 and the NACRE rate, as a function of tempera- 
ture. Hatched areas: values of T c where H and He burning phases 
take place in a 15 Mq model at Z = 0. 

The variation of the partial widths of both reactions have 
been computed in Appendix B and are depicted in Fig. IB. II 
Together with the results of the previous section and the de- 
tails of the Appendix B, we can compute the 3a-reaction rate 
as a function of temperature and Snn- This is summarized in 
Fig. [3] which c ompares the rate fo r different values of 5nn to the 
NACRE rate dAngulo et al.l fl999). which is our reference when 
no variation of constants is assumed (i.e. 5nn - 0). One can also 
refer to Fig. IB. 21 which compares the full numerical integration 
to the analytical estimation which turns out to be excellent in 
the range of temperatures of interest. As one can see, for positive 
values of 6nn, the resonance energies are lower, so that the 3a 
process is more efficient (see Appendix IB1. 

Let us compare the result of Fig. [3] which gives y = 
log[/l3 Q ,(^iv)//l3 U ,(0)] to a simple estimate. Using the analytic 
expression (f2]i for the reaction rate, valid only for a sharp res- 
onance, y is simply given by 



1 



y = 



In 10 



(15) 



where the sensitivity Sg m = dln/l3 Q ./dln5iVA' is given, from 
Eq. Q3), by s Srm = 6 NN x (32.62MeV)/Jtr. We conclude that 



y = 1.644 xfe)(-J-) 

y \io- 3 mo 8 K/ 



(16) 



This gives the correct order of magnitude of the curves depicted 
in Fig.[3]as well as their scalings with 6nn and with temperature, 
as long as Tg > 0. 1 . At lower temperatures differences arise from 
the fact that the analytical expression for the reaction rate is no 
longer accurate (see Appendix IB1. 

The sensitivity to a variation of the intensity of the N-N in- 
teraction arises from the fact that dQ aaa /d6NN ~ ^0 2 Q„ a „. The 
fact that the typical correction to the resonant energies is of or- 
der 10 MeV (x 6nn) compared to the resonant energies them- 
selves which are of order 0. 1 MeV, allows one to put relatively 



strong constraints on any variation. This is reminiscent of the 
case of the resonance producing an excited state of 150 Sm of im- 
portance in setting c onstraints on the variation of couplings us- 
ing the Oklo rea c tor IShlvakhterlll976l;lDamour & Dysonf l 996; 
iFuiii etalj|2000l lOlive et alj 120021; iPetrov et all 120061) . In that 
case, the resonant energy is 0. 1 eV compared to corrections of 
order 1 MeV due to changes in the fine structure constant, lead- 
ing to limits on Aa em /a em of order 10~ 7 . 

3.4. Using the Deuterium binding energy as a link to 
fundamental constants 

The nuclear model described above introduces the parameter 
Snn which is itself not directly related to a set of fundamental 
constants such as gauge and Yukawa couplings. In order to make 
such a connection, we make use of previous analyses relating the 
deuterium binding energy Bp to fundamental constants. 

Using a potential model, the dependence of Bp 
on the nucleon, cr-meson and tii-meson has been esti- 
mated 



2003 



(Flambaum & Shurvak 
Flambaum & Shurvak 



Coc et al 



2002; Dmitriev & Flambaun 
2003; Dmitriev et al.1 2004 



120071: iDamour & Donoghud 120081) . Furthermore, 
using the quark matrix elements for the nucleon, variations 
in Bo can be related to variations in the light quark masses 
(particularly the strange quark) and thus to the corresponding 
quark Yukawa couplings and Higgs vev, v. The remaining 
sensitivity of B p to a dim e nsion ful quantity is ascribed to the 
QCD scale A. In lCoc etal] d2007l) . it was concluded that 



B n A I v h< 



(17) 



Eq. ([8j can then link any constraint on 6nn to the three funda- 
mental constants (h s , v, A). 

Further relations are possible in the context of unified theo- 
ries of gauge interactions. From the low energy expression for 
Aqcd, 



A 



m c nib nit \ 



2/27 



exp - 



In 



9a s (M) 



(18) 



one can determine the relation between the changes in A and 
the g au ge couplings an d quark masses (ICampbell & Olive] 
1995: iLangackeretall 120021: iDent & Fairbairnl [2003; 
ICalmet & F ritzsch 2002; lDamour et al.ll2002al) . 



AA 



= *^L + ! 3^ + ^ + ^ + ^|. (19) 
A a em 27 \ v h c hb h, j 

Typical values for R are of order 30 in many grand unified the- 
ories, but there is cons iderable model dependence in this coeffi- 
cient dDine et al]|2003l) . 

Furthermore, in theories in which the electroweak scale is 
derived by dimensional transmutation, changes in the Yukawa 
couplings (particularly the top Yukawa) leads to exponentially 
large changes in the Higgs vev. In such theories, 



Av Ah 
— ~S — 
v h 



(20) 



with S ~ 160, though there is considerable model dependence 
in this value as well. 

Finally, using the relations in Eqs. dT9b and d20l i. we can 
write 



AB D Ah Aa em 

— - = -13(1+5) — + 18/? — — 

B D h a em 



(21) 
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If in addition, we relate the gauge and Yukawa couplings through 
Ah/h = (l/2)Aa em /a em , we can further write, 



AB L 

Bn 



= -[6.5(1 + S)- 18/?]- 



Aa» 



(22) 



An alterna t ive i n vestigation ( Yoo & Scherrerl [2003 ; 
Epelbaum et al. I 120031: 1 ean & Savage! 120031) suggests a 
large dependence of Bo on the pion mass m„, 



AB L 
Bn 



Am,, 



(23) 



where r is expected to range between 6 and 10. Again, this al- 
lows one to related Bp, and thus 6nn to (h, v, a em ). 

As these two examples demonstrate, the main problem arises 
from the difficulty to determine the role of the QCD parameter 
in low energy nuclear physics. They show, however, that such a 
link can be drawn, even though it is strongly model-dependent. 

4. Stellar implications 

The Geneva stellar code was adapted to take into account the re- 
action rates compute d above. The version of the code we use is 
the one described in Ekstror rTet al] (120081) . Here, we only con- 
sider models of 15 Mq and 60 Mq without rotation and as- 
sume an initial chemical composition given by X = 0.7514, Y = 
0.2486 and Z = 0. This corresponds to the BB N abundance of 
He at the baryon density determined by WMAP (Komatsu et al. 
2009) and at zero metallicity as is expected to be appropriate for 
Population III stars. For 16 values of the free parameter 6nn m 
the range -0.009 < 6nn ^ +0.006, we computed a stellar model 
which was followed up to the end of core He burning (CHeB). 
As we will see, beyond this range in 5nn, stellar nucleosynthesis 
is unacceptably altered. Note that for some of the most extreme 
cases, the set of nuclear reactions now implemented in the code 
should probably be adapted for a computation of the advanced 
evolutionary phases. 

Focusing on the limited range in 5nn will allow us to study 
the impact of a change of the fundamental constants on the pro- 
duction of carbon and oxygen in Pop III massive stars. In this 
context, we recall that the observations of the most iron poor 
stars in the halo offer a wonderful tool to probe the nucleosyn- 
thetic impact of the first massive stars in the Universe. Indeed 
these halo stars are believed to form from material enriched by 
the ejecta of the first stellar generations in the Universe. Their 
surface chemical composition (at least on the Main Sequence), 
still bear the mark of the chemical composition of the cloud from 
which they formed and thus allow us to probe the nucleosyn- 
thetic signature of the first stellar generations. Any variation of 
the fundamental constants which for instance would prevent the 
synthesis of carbon and/or oxygen would be very hard to con- 
ciliate with present day observations of th e most iron poor stars . 
For instance the tw o most iron poor stars (Christ lieb et al.f 2004; 
Frebel et al. 2008) both show strong overabundances of carbon 
and oxygen with respect to iron. 

Our results for 15 Mq and 60 Mq stars are presented in § 14.11 
and § |4.2| respectiveiv. 

4.1. 15 Mq mass star 

Figure |4](fe/f) shows the HR diagram for the models with 6nn 
between -0.009 et +0.006 in increments of 0.001 (from left to 
right). Once the CNO cycle has been triggered (see below), the 
Main Sequence (MS) tracks are shifted towards cooler T e g for 



increasing 6nn- There is a difference of about 0.20 dex between 
the two extreme models. Figure \4\(right) shows the central tem- 
perature at the moment of the CNO-cycle ignition (lower curve). 
On the zero-age main sequence (ZAMS), the standard (6nn = 0) 
model has not yet produced enough 12 C to be able to rely on the 
CNO cycle, so it starts by continuing its initial contraction until 
the CNO cycle ignites. In this model, CNO ignition occurs when 
the central H mass fraction reaches 0.724, i.e. when less than 
3% of the initial H has been burned. Models with 6nn < (i.e. 
a lower 3a rate) yield a phase of contraction which is longer for 
lower 5nn (i.e. larger \5nn\Y m these models, the less efficient 
3a rates need a higher T c to produce enough 12 C for triggering 
the CNO cycle. Models with 6nn > (i.e. a higher 3a rate) are 
directly sustained by the CNO cycle on the ZAMS: the star can 
more easily counteract its own gravity and the initial contraction 
is stopped earlier, so H burning occurs at lower T c and p c (Fig. [4] 
right), i.e. at a slower pace. The MS lifetime, tms, is sensitive to 
the pace at which H is burned, so it increases with 6nn- The rel- 
ative difference between the standard model MS lifetime tms at 
6 NN = and t M s at 6 NN = -0.009 (+0.006) amounts to -17% 
(+19%). 

While the differences in the 3a rates do not lead to strong ef- 
fects in the evolution characteristics on the MS, the CHeB phase 
amplifies the differences between the models. The upper curve 
of Fig. ^(right) shows the central temperature at the beginning 
of CHeB. There is a factor of 2.8 in temperature between the 
models with 6nn - -0.009 and +0.006. To get an idea of what 
this difference represents, we can re late these temperat ures to the 
grid of Pop III models computed bv lMarigo et al .1 (12001 1) . The 15 
Mq model with 6nn = -0.009 starts its CHeB at a higher tem- 
perature than a standard 100 Mq of the same stage. In contrast, 
the model with 6f/N = +0.006 starts its CHeB phase with a lower 
temperature than a standard 12 Mq star at CNO ignition. Table|2] 
presents the characteristics of the models for each value of 6nn 
at the end of CHeB. From these characteristics, we distinguish 
four different cases (see the last column of Tableland Fig . [5j : 



I In the standard model and when 6mn is ver y close to 0, 12 C 
is produced during He burning until the central temperature 
is high enough for the 12 C(a,y) I6 reaction to become effi- 
cient: during the last part of the CHeB phase, the I2 C is pro- 
cessed into 16 0. The star ends its CHeB phase with a core 
composed of a mixture of 12 C and 16 (see the top left panel 
of Fig. 0. 

II If the 3^ rate is weakened (-0.005 < 6 NN < -0.002), 12 C is 
produced at a slower pace, and T c is high from the beginning 
of the CHeB phase, so the 12 C(a,y) I6 reaction becomes 
efficient very early: as soon as some 12 C is produced, it is 
immediately transformed into 16 0. The star ends its CHeB 
phase with a core composed mainly of 16 0, without any 12 C 
and with an increasing fraction of 24 Mg for decreasing 6nn 
(see the bottom left panel of Fig. |5). 

III For still weaker 3a rates (6nn ^ -0.006), the central temper- 
ature during CHeB is such that the 16 O(ar,y) 20 Ne(Q',y) 24 Mg 
chain becomes efficient, reducing the final 16 abundance. 
The star ends its CHeB phase with a core composed of nearly 
pure 24 Mg (see the bottomright panel of Fig.|5]l. Because the 
abundances of both carbon and oxygen are completely neg- 
ligible, we do not list the irrelevant value of C/O for these 
cases. 

IV If the 3a rate is strong (6 NN > +0.003), 12 C is very rapidly 
produced, but T c is so low that the I2 C(a,y) I6 reaction can 
hardly enter into play: 12 C is not transformed into lfj O The 
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Fig. 4. Left panel: HR diagrams for 15 Mq models with 6nn - (thick black), and ranges from left to right from -0.009 to +0.006 
in steps of 0.001 (using same colour code as Fig. [3). Right panel: The central temperature at CNO ignition (circles) and at the 
beginning of CHeB (squares) as a function of 5nn- The labels on the CNO-ignition curve shows the central H mass fraction at that 
moment. Note that above 5nn — 1-0.001, CNO ignition occurs on or before the ZAMS. 



star ends its CHeB phase with a core almost purely com- 
posed of 12 C (see the top right panel of Fig. [5J. 

These results are summarized in Fig.|6]which shows the compo- 
sition of the core at the end of the CHeB phase. One can clearly 
see the dramatic change in the core composition as a function 
of 6nn showing a nearly pure Mg core at large and negative 
6nn, a dominantly O core at low but negative 6nn, an d a nearly 
pure C core at large and positive 6nn- These results are quali- 
tatively consistent with those found by Sch lattl et al.l (120041) for 
Population I type stars. Note that their cases with A£r = +100 
keV correspond roughly to our 6nn ~ +0.005. 

Table [2] s hows also the core size at the end of CHeB. As in 
iHeger et al l (120001) the mass, Mqo, is determined as the mass 
coordinate where the mass fraction of 4 He drops below 10" 3 . 
The mass of the CO core increases with decreasing Snn, the in- 
crease amounting to 8% between 6nn = +0.006 and -0.009. This 
effect is due to the higher central temperature and greater com- 
pactness at low 6nn . The same effec t was found by other authors 
dSchlattl et al.l2004tlTur et al.l2007h . As shown by these authors, 
this effect is expected to have an impact on the remnant mass and 
thus on the strength of the final explosion. 

4.2. 60Mq mass star 

As it is widely believed that Pop III stars are massive, we next 
present results for 60 Mq models (at Z = 0). The characteris- 
tics of these models for different values of 6nn are collected in 
Table [3] 

Figure [7] shows the HR diagram for our 60 Mq models. 
During the MS, the shift of the tracks in T e g are slightly reduced 
compared to the 15 Mq models: by 0. 18 dex. Also, all the 60 Mq 
models are directly sustained by the CNO cycle on the ZAMS, so 
the tracks are just shifted regularly, without affecting the shape 
of the tracks. During CHeB, however, the behavior we described 



Table 2. Characteristics of the 15 Mq models with 6nn ran gi n g 
from -0.009 to +0.006 at the end of core He burning. The MS life- 
time, tms, and the core He-burning duration, TcHeB, are expressed 
in Myr, the CO-core mass, Mco, is in Mq, X{C) is the central 
value for the carbon mass fraction and C/O is the ratio of the car- 
bon to oxygen mass fractions. 
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C/O 


case 


-0.009 


8.224 


1.344 


3.84 


4.4e-10 




III 


-0.008 


8.285 


1.276 


3.83 


2.9e-10 






-0.007 


8.308 


1.200 


3.38 


8.5e-10 






-0.006 


8.401 


1.168 


3.61 


4.2e-07 






-0.005 


8.480 


1.130 


3.59 


5.9e-06 


3.0e-05 


II 


-0.004 


8.672 


0.933 


3.60 


3.2e-05 


5.2e-05 




-0.003 


8.790 


0.905 


3.60 


1.3e-04 


1.7e-04 




-0.002 


9.046 


0.892 


3.61 


5.6e-04 


6.4e-04 




-0.001 


9.196 


0.888 


3.70 


0.013 


0.014 


I 





9.640 


0.802 


3.65 


0.355 


0.550 




+0.001 


9.937 


0.720 


3.61 


0.695 


2.278 




+0.002 


10.312 


0.684 


3.62 


0.877 


7.112 




+0.003 


10.677 


0.664 


3.62 


0.958 


22.57 


IV 


+0.004 


10.981 


0.659 


3.62 


0.981 


52.43 




+0.005 


11.241 


0.660 


3.61 


0.992 


123.9 




+0.006 


11.447 


0.661 


3.55 


0.996 


270.2 





" mass coordinate where the abundance of 4 He drops below 10 
b Here, X(C) = X( [2 C) 



for the 15 Mq models with 6nn < is more pronounced in the 
case of the 60 Mq models: 12 C and ls O are already exhausted 
at the end of CHeB (case III) for 6nn < -0.005. This can be 
understood because the 12 C(a,y) 16 0, the 16 O(a,y) 20 Ne and the 
20 Ne(a,y) 24 Mg reaction rates, are a factor of 10 to 100 higher 
than the 3a rate when log T c * 8.48, i.e. when there is still about 
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Fig. 5. The evolution of the central mass fraction for the main 
chemical species inside the core of the 15 Mq models during 
core He burning. Top left: standard case (representative of case 
I, see text), bottom left: 5nn = -0.003 (case II), bottom right: 
S m = -0.007 (case III), top right: 6 NN = +0.005 (case IV). 
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Fig. 6. The composition of the core at the end of the central He 
burning in the 15 Mq models as a function of Snn- 



5% of helium in the core. Instead of a CO core, these models are 
left with an almost pure 24 Mg core. 

For the 60 Mq models with S NN > 0, there is still a reason- 
able abundance of oxygen up to 5nn = +0.003. At higher values 
of 6nn, we are again left with a nearly pure carbon core. The 
model with 6^ = +0.006 has proven to be very difficult to fol- 
low at the end of CHeB and was stopped before complete He 
exhaustion. The results for the 60 M0 models are summarized 
in Fig. [8] which shows the composition of the core at the end of 
the CHeB phase. As in the case of the 15 M0 models, one can 
clearly see the strong dependence of the core composition on 
•Saw- 



Table 3. Characteristics of the 60 Mq models with 5nn ranging 
from -0.009 to +0.006. tms and TcHeB are expressed in Myr, Mco 
is in Mq, X(C) is in mass fraction. The mark "75" after Mco 
indicates that it has been calculated as the mass coordinate where 
the C+O abundance rises above 0.75. The mark "75Mg" after 
Mco indicates that it is the 24 Mg abundance which has become 
higher than 0.75. 
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II 


-0.003 
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2.1e-04 
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3.338 


0.328 
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9.3e-04 


1.4e-03 




-0.001 


3.379 
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18.40' 5 


0.008 


0.009 
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3.418 


0.299 


21.37 


0.163 


0.200 




+0.001 


3.458 


0.267 


21.59 


0.513 


1.062 




+0.002 


3.495 


0.244 


21.14 


0.761 


3.193 




+0.003 


3.534 


0.259 


21.14 


0.899 


9.233 




+0.004 


3.571 


0.265 


19.86' 3 


0.952 


19.89 


IV 


+0.005 


3.607 


0.808 


27.22 


0.963 


32.85 




+0.006 


3.644 













The effect of varying 6mn on the core size is less clear in 
the case of the 60 Mq models. In some cases, the model under- 
goes a CNO boost in the H-burning shell during CHeB, which 
reduces the core masfl The occurance of the boost does not fol- 
low a clear trend with 6nn- It appears on the HR diagram as a 
sudden drop in luminosity and effective temperature in the red- 
wards evolution during CHeB (see Fig. [4] left). 



4.3. Limits on the variation of the fundamental constants 

All of the models considered were followed without any numer- 
ical or evolutionary problem through the MS. The differences 
in lifetimes and tracks during this phase are not constraining 
enough to allow the exclusion of some range in 6nn between 
-0.009 and +0.006. However, the CHeB phase amplifies these 
differences. 

At the end of CHeB, the models with 6 NN < -0.005 for the 



15 Mq model and 6, 



• m < -0.004 for the 60 M Q model have vir- 
tually no 12 C in the core, which means that the "standard" suc- 
cession of stellar evolution burning phases will not be respected 
(see the bottom right panel of Fig. [5). These models are also de- 
void of I6 or 20 Ne as well, leaving us with a nearly pure 24 Mg 
core. Note that at this phase, the central temperature is close to 
that which would allow the 24 Mg(y,o-) 20 Ne or 24 Mg(a,y) 28 Si re- 
actions to take place. Therefore, there is a possibility that the 
nucleosynthetic chain could go on despite its strange evolution. 
However, the Geneva code is developed to follow the standard 
phases of stellar evolution, making it necessary to be modified 
before being able to follow further the evolution of these odd 
objects. 

The models with 6nn between -0.002 and -0.005 (between 
-0.002 and -0.004 for the 60 Mq model) end the CHeB phase 
with a central abundance of I2 C between 10 and 10 -7 , which 
means that the central C-burning phase will be extremely short. 
The 2() Ne abundance at that stage is comprised between 0.04 and 



2 W e refer the interested reader to Hirschi (20071) or lEkstrom et al.l 
( 2008) for a more detailed description of this phenomenon. 
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Fig. 7. Left panel: HR diagrams for 60 Mq models with Snn = (thick black), and ranges from left to right from -0.009 to +0.006 
in steps of 0.001 (using the same colour code as in Figs.[3]and @}. Right panel: The central temperature at CNO ignition (circles) and 
at the beginning of core He burning (squares) as a function of Snn- In the 60 M0 models, CNO-ignition occurs while the hydrogen 
mass fraction at the centre is still equal to its initial value (0.75). 
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Fig. 8. The composition of the core at the end of the central He 
burning in the 60 Mq models as a function of Snn- 



0.10, so there will be a short phase of neon photodisintegration. 
Moreover, the 16 abundance ranges between 0.94 and 0.44 so 
the oxygen fusion phase will be almost normal. While the suc- 
cession of the burning phases seems preserved, one can however 
suppose that these models will present very different yields than 
the standard case with Snn = 0. This point could be the subject 
of a future study. It is interesting to note here that since the C- 
burning phase is very short (due to the very low 12 C abundance 
at the end of CHeB), the model will not have much time to lose 
entropy by neutrinos losses. We can suppose that the iron core 
will be hotter and bigger, so the remnant could be a black hole 



instea d of a neutron star ( Wooslev & Weaver! 1986 : Schlatt let al.l 
12004 . 

The models with Snn > end the CHeB phase with larger 
and larger 12 C abundances for increasing Snn- The carbon burn- 
ing phase will thus be much longer for these models which will 
lose a lot of energy through neutrino emission. A more suspi- 
cious feature is that the ls O production becomes negligible or 
even null for Snn ^ +0.003 (see the top right panel of Fig. [5]) 
(> +0.004 for the 60M Q model). Normally the bulk of the 
16 production occurs during CHeB: during C burning, the lfi O 
abundance is reduced by I6 O(a,y) 20 Ne, and during Ne burning, 
only a small fraction is produced by the photo-disintegration re- 
action 20 Ne(y,ar) I6 O. It would thus mean that such stars do not 
produce any I6 0. This would pose difficulties for explaining the 
high O overabundances obser ved in extremely iron-poor stars 
found in the Galactic halo (see lFrebel et al. 2008). 

From the preceding discussion, if we exclude a core com- 
posed exclusively of 24 Mg (case III), we must reject Snn < 
-0.005 for the 15 Mq model. If we consider that a core only 
composed of 12 C is not acceptable either (case IV), we must re- 
ject Snn > +0.002. If we consider that a reasonable value of 
C/O must lay close to unity, we must also reject case II and the 
allowed range for Snn is further restricted to -0.001 to +0.002. 
Similarly for the 60 Mq model, excluding cases III and IV leads 
to a limit -0.004 < Snn ^ +0.003. The more stringent condition 
on C/O ~ 1 leads to -0.001 < S NN < +0.003. 

5. Discussion 

As we have seen in the previous sections, the extreme sensitivity 
of the triple a process to the resonant energy of the Hoyle state 
can lead to very different histories for massive Population III 
stars. In particular, we have shown that very slight variations in 
the nucleon-nucleon potential (of order a few xl0~ 3 ) can lead 
to very different core compositions at the end of CHeB. We 
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identified two cases (III and IV) corresponding to nearly pure 
24 Mg or pure 12 C cores. These cases were present in both the 
15 and 60 M0 models studied. Below 6nn = -0.005, the stars 
end the CHeB phase with a core that is almost completely de- 
prived of carbon, oxygen and neon. This is due to the 12 C pro- 
duction by the 3a reaction becoming extremely weak compared 
to the 12 C(a,y) 16 reaction (for which we have used the rates of 
iKunz et alJl2002h . As soon as a little amount of 12 C is produced, 
it is transformed into 16 0, which in turn is transformed into 20 Ne 
and then 24 Mg due to the high temperature and density at which 
He burning occurs in these models. Above Snn = +0.002, the 
models end the CHeB phase devoid of ls O. 

We have checked the limiting values for 6nn variation with 
stellar models in two different mass domains. For the 15 Mq 
models, the lower limit is slightly larger than for the 60 Mq 
models. This is due to the fusion phases occurring at higher T c 
in the more massive stars, at conditions where the 12 C(ff,y) 16 0, 
l6 O(a,y) 20 Ne and 20 Ne(o',y) 24 Mg reaction rates are largely dom- 
inant over the 3a rate. A weak 3a reaction is a bigger handicap 
in the high mass domain. In contrast, the upper limit is larger for 
the 60 Mq models, because the 3a reaction rate is a little less 
extreme at higher T c . 

Excluding these cases allows us to set a relatively conserva- 
tive limit on 6nn> 



0.004 < 5 NN < +0.002 . 



(24) 



A more aggressive limit would also exclude case II in which 
CHeB ends with a 16 and 20 Ne core with little or no 12 C. In this 
case, one could argue 



- 0.001 < 6 NN < +0.002 . 



(25) 



For the remainder of the discussion, we will restrict our attention 
to the weak limit ( l24b . as our conclusions can be easily scaled to 
the stronger limit. 

The limit in Eq. ( l24l i stems directly from the variation in 
Qaaa- Excluding regions III and IV amount to limiting Q aaa to a 
range 0.3142 - 0.5100 Me V, or 



017< A£W <Q 34 

Qaaa 



(26) 



As discussed in 33.21 a variation in 6nn will result in a vari- 
ation in the deuterium binding energy. Using Eq. (0, the bound 
d24l i thus becomes 



0.023 < < +0.011 . 
Bn 



(27) 



In principle, one would like to next convert the limit on 6nn or 
B D into a limit on the fundamental constants. Unfortunately, as 
we have argued earlier, i) the direct limit from the triple a pro- 
cess based on 6 a is far weaker than that due to 6nn and ii) in the 
absence of some guiding theory of unification, we can not relate 
the variation of Ba directly to a var iation in a^m- How ever, as dis- 
cussed above and in more detail in lCoc et al. (2007), we can use 
gauge coupling unification to relate a variation in Bo to a varia- 
tion in a em through A. Ignoring first any variation in the Yukawa 
couplings and Higgs vev, thus using ABq/Bo = l8RAa em /a em , 
with R — 36 as is expected in the simplest grand unified theories, 
we obtain 



3 .5xl0- 5 <^<+1.8xl0- 5 . 



(28) 



If we further assume the relations between gauge and Yukawa 
couplings and use Eq. d22l . the limit, though more speculative, 



is actually weakened by a factor of about 2 due to the partial 
cancellation between the gauge and Yukawa contributions to B D . 

The limits on the variation of the fine structure constant de- 
rived above corresponds to a variation in a em between the present 
time and a period around a redshift z ~ 15-20 where the 
Population III stars would have been present. These values are 
compatible with the similar limits (also assuming gauge cou- 
pling unification) on the variation of a em at a redshift of 10 10 
from BBN predictions. They are larger by a factor of 10 than 
the values found in the claimed detections dWebbetal]|200U 
iMurphv et ail I2003L 120071) or a factor of 10 weaker than the 
limits from the non-detection (IChand et alj|2004t ISrianand et al.l 
2004; lOuast et al.ll2004t ISrianand et alj|2007l) of a variation in 
a em from quasar absorption systems at redshifts z < 3.5. 

We remind the reader, that in the present work, the varia- 
tion in Snn is only taken into account for the 3a reaction. If the 
other rates are also affected, the limits found here would poten- 
tially have to be revised, because they have been determined by 
anomalies in the evolution that are due to a competition between 
the efficiency of the various rates. However, being a resonant 
reaction, t he 3a reaction is expecte d to be the most sensitive. 
Following bberhummer et all (l2000h . the 16 O(ar,y) 20 Ne reaction 
is not expected to be sensitive to variations in a em , while the 
l2 C(a,y) lb O could be more affected by such variations because 
of subthresholds in the 16 nucleus. According to the same au- 
thors, this last reaction is expected to be strengthened by a weak- 
ening of the nucleon-nucleon interaction. In this case, the effects 
described for cases II and III would be more dramatic than the 
ones presented here, and so the limits might be tighter. 

We conclude by asking: is it reasonable to exclude values of 
5nn using nucleosynthetic constraints from stellar models? The 
criteria that we applied assumes the possibility for a "normal" 
succession of burning phases (H — > He — » C — > Ne — » O — » Si). 
Although we have not done so here, a modified code including 
"non standard" fusion phases, would allow us to follow these 
models further. It is expected that the resulting yields would 
present large anomalies. Given the current state of abundance 
determination in extremely metal poor stars, it is highly improb- 
able that the first stars would not produce fair amounts of 12 C 
and lfi O. The conservative case seem thus to offer a reasonable 
limit on the variations of the fundamental constants. 



Appendix A: Details on the Microscopic model 

Here, we provide some technical details about the micro- 
scopic calculation used to determine the 8 Be and 12 C bind- 
ing energies. This calculation is based on the description of 
the nucleon-nucleon in teraction by the Minnesota (MN) force 
dThompson et al.|[T977l) . adapted to low-mass systems. 

The nuclear part of the interaction potential Vn between nu- 
cleons i and j is given by 



V NiJ (r) 



where r ■ 



V R (r) + 1(1 + P%)V t (r) + 1(1 - Pfj)V s (r) 



2 
1 

2" ' 2 



1 1 

u + -(2 - u)P r ij 



(Al) 



-r ,| and Pf. and P r ., are the spin and space exchange 

J 'j './ 

operators, respectively. The radial potentials Vfi(r), V s (r), V t (r) 
are expressed as Gaussians and have been optimized to repro- 
duce various properties of the nucleon-nucleon system, such as 
the deuteron binding energy at 6nn — 0, or the low-energy phase 
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shifts. They have been fit as (Thompson et all! 9771) 

V R (r) = 200 exp(-1.487r 2 ) 
V s (r) = -91.85 exp(-0.465r 2 ) 
V,(r) = -178 exp(-0.639r 2 ) 



(A2) 



where energies are expressed in MeV and lengths in fm. 

In Eq. dAU . the exchange-admixture parameter u takes stan- 
dard value m = 1 , but can be slightly modified to reproduce im- 
portant properties of the A-nucleon system (for example, the en- 
ergy of a resonance). This does not affect the physical proper- 
ties of the interaction. The MN force is an effective interaction, 
adapted to cluster models. It is not aimed at perfectly repro- 
ducing all nucleon- nucleon properties, as realistic forces used 
in ab initio models (iNavratil et al.ll2009h . where the cluster ap- 
proximation is not employed. The potentials are expressed as 
Gaussian factors, well adapted to cluster models, where the nu- 
cleon orbitals are also Gaussians (IWildermuth & Tangill977h . 

The wave functions (|6]l are written in the Resonating Group 
Method (RGM) which clearly shows the factorization of the sys- 
tem wave function in terms of individual cluster wave func- 
tions. In practice the radial wave functions are expanded over 
Gaussians, which provides the Gen erator Coordinate Method 
(GCM), fully equivalent to the RGM (IWildermuth & Tangll977h 
but better adapted to numerical calculations. Some details are 
given here for the simpler two cluster case. The radial function 
g J 2 Mn (p) is written as a sum over Gaussian functions centered at 
different values of the Generator Coordinate R„. This allows us 
to write the 8 Be wave function (0 as 



xuJMn 



(A3) 



where <S) JMn (R n ) is a projected Slater determinant. 

This development corresponds to a standard expansion on 
a variational basis. The binding energies E Jn of the system are 
obtained by diagonalization of 



2 [H J \R n ,R w ) - E J "N J "(R n ,R„,)]f J \R„) = 0, 



(A4) 



the time scale for alpha capture on 8 Be is negligible compared 
to its lifetime against alpha decay. The rate is calculated as in 
NACRE for the resonance of interest: 



N 2 A (o-v) aaa 



3N A 



Ha 



2nk B T 



3/2 



f 
Jo 



^lexp(-E/k B T)N A (o-vf^ 
T a (E) 



EdE, 



(Bl) 



where fj, aa is the reduced mass of the a + a system, and E is the 
energy with respect to the a + a threshold. The elastic cross sec- 
tion of a + a scattering is given by a Breit-Wigner expression: 



( w\ - 71 r 2 (£) 

0- aa{ E) - k 2" (E _ ER)) 2 +r 2 (E)/4 > 



(B2) 



where k is the wave number, E R =E R ( & Be), Y a =Y a (*~Be), wis a 
statistical factor (here equal to 2 to account for identical particles 
with spin zero). 

The N A {(Tv) a Be rate assu mes that 8 Be has been formed at an 
energy E different from E» Be (iLanganke et al.|[T986l) . This rate is 
given by 



N A (crvy 



a'Be 



= N A 



8tt / /v Be 



ju' \2nk B T 



3/2 



f 

Jo 



cr a8Be (£"; E) exp(-E'/k B T) E' dE' , 



(B3) 



where /^ a s Be is the reduced mass of the a + 8 Be system, and 
E' is the energy with respect to its t hreshold (whic h varie s 
with the formation e nergy E). As in Nomo toet al.l (Il985l) : 
Langa nke et al.l d 19861) . we parametrize cr a s Be (E'; E) as 



Be (E';E) 



nn 2 



2fi a s Be E' 



where the overlap and hamiltonian kernels are defined as 



N jK (R n ,R n> ) 
H Jn (R„,R„,) 



(® Jn (R„)\<l> J *(R n ,)), 
(® J,t (R n )\H\® J *(R n ,)). 



(A5) 



The Hamiltonian H is given by Eq. ([3]). Standard techniqu es ex- 
ist for the evaluation of these many -body matrix elements (Brink 
Il966t) . The choice of the nucleon-nucleon interaction directly af- 
fects the calculation of the hamiltonian kernel, and therefore of 
the eigenenergy E Jn . 

For three-body wave function, the theoretical develop- 
ments are identical, but the presentation is more complicated 
due to the presence of two relative coordinates (p,R). The 
problem is addressed by using th e hyperspherical formalism 
dKorennov & Descouvemontl20 04 ) . 

Appendix B: Reaction rates and numerical 
integration 

To take into account the (energy dependent) finite widths of the 
two resonances involved in this two step process, one has to per- 
form numerical integ rati ons as was done in NA CRE following 
iNomoto et ail dl985f) and lLanganke etail dl986l) . Here, the con- 
dition of thermal equilibrium is relaxed, but it is assumed that 



r a (E')T y (E' + E) 



[E' - £ S ( 12 C) + E - £ s ( 8 B e )] 2 + \T(E>; E) 2 



(B4) 



where the partial widths are those of the Hoyle state and in par- 
ticular, F = F a ( l2 C) + F r ( 12 C). The various integrals are calcu- 
lated numerically. The experimental widths at resonance energy 
can be found in Table Q] 

However, one must include i) the energy dependence of those 
widths, away from the resonance energy and ii) the variation of 
the widths at the resonant energy when this energy changes due 
to a change in the nuclear interaction. 

The energy dependence of the particle widths r o (E) is given 

by: 



T a (E) = T a (E R ) 



P((E,R C ) 
P((E R ,R c y 



(B5) 



where Pc is the penetration factor associated with the relative 
angular momentum t (0 here) and the channel radius, R^ The 
penetration factor is related to the Coulomb functions by: 



Pr(E,R) = 



Fj(r,,p) + Oj(r},p) 



(B6) 



3 We choose R c = 1.3 (A + XL ) fm, for nuclei A i and A 
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Fig. B.l. The partial widths of 8 Be and I2 C as a function of 6 



DC 

° -5 
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Fig. B.2. The 4 He(aa, y) 12 C reaction rate as a function of tem- 
perature for different values of dw Solid (dashed) lines repre- 
sent the result of the numerical calculation (analytical approxi- 
mation) with 6nn-0 (black), 6nn > (red) and 6nn < (blue). 
(For the negative value of S^n, the larger difference is caused by 
the failure of the numerical integration and the analytical solu- 
tion is preferred (see text).) 



where p - kR and 



V = 



ZiZ 2 a e , 



v/c 



(B7) 



is the Sommerfeld parameter. 

For radiative capture reactions, the energy dependence of the 
gamma width T y {E) is given by: 



r y (E)<xa em E 



(B8) 



where A is the multipolarity (here 2 for E2) of the electromag- 
netic transition. 

The relevant widths as a function of 6nn are given in 
Figure IB. 11 They are directly linked to the resulting change of 
£ R ( 8 Be) and£ s ( 12 C). 

The radiative width, r y ( 12 C) with its E 5 energy dependence 
shows little evolution. (The energy of the final state at 4.44 MeV 
is assumed to be constant). In contrast, the 8 Be alpha width un- 
dergoes large variations due to the effect of Coulomb barrier 
penetrability. Note that compared to these variations, those in- 
duced by a change of a em in the Coulomb barrier penetrability 
(Eqs. |B7| |B61 > and T y are considerable smaller. 

Numerical integration is necessary at low temperature as 
the reaction takes place through the low energy wing of reso- 
nances. It takes even more relative importance, at a given tem- 
perature, when the resonance energy is shifted upwards. On 
the other hand, when 6nn increases, the resonance energies de- 
crease, and the r„( 8 Be) becomes so small that the numerical in- 
tegration becomes useless and soon gives erroneous results be- 
cause of the finite numerical resolution. For this reason, when 
r„.( 8 Be) < 10~ 8 MeV, we use instead the Saha equation for the 
first step and the sharp resonance approximation for the second 



step, i.e. Eq. © when r a ( 12 C) < 10" 8 MeV. (Note that for high 
values of 6nn, the condition T y <k T a does not hold anymore and 

r * r r .) 

At temperatures in excess of 2, one must include the 
contribution of the higher 12 C levels like the one observed by 
Fvnb o et all d2005). As this is not of importance for this study, 
we just added the contribution given by the last terms in the 
NACRE analytical approximation and neglected any induced 
variation. 

Figure lB~2l shows the numerically integrated 4 He(ao', y) 12 C 
reaction rates for different values of 6nn compared with the an- 
alytical approximation (Eq. (|2]i). The difference is important at 
low temperature and small 6mn values but becomes negligible 
for &nn k. 0. At the highest values of 6nn we consider, the nu- 
merical calculation uses the Saha equation for the first step but 
the total widths of the I2 C level becomes also too small to be 
accurately numerically calculated: we use Eq. instead. 

The 8 Be lifetime w.r.t. alpha decay, (h/r a ( & Be)), exhibits 
the opposite behavior indicating that for large values of 6nn 
it becomes stable. Before that, its lifetime is so long that the 
4 He(ao-, y) 12 C reaction should be considered as a real two step 
process with 8 Be included in the network as the assumption that 
alpha decay is much faster than alpha capture may not hold any- 
more. Fortunately, our network calculations shows that this sit- 
uation is encountered only for 6nn S 0.006 for the temperatures 
and densities considered in our stellar evolution studies. 



Acknowledgements. This work was partly supported by i) CNRS PEPS/PTI, 
ii) CNRS PICS France-USA and by DOE grant DE-FG02-94ER-40823 at the 
University of Minnesota. 



14 



Ekstrom et al.: Effect of the variation of fundamental constants on Pop III star evolution 



References 

Adams, F. C. 2008, JCAP, 0808, 010 
Ajzenberg-Selove, F. 1990, Nucl. Phys. A, 506, 1 
Angulo, C, Arnould, M., Rayet, M., et al. 1999, Nucl. Phys. A, 656, 3 
Ashenfelter, T., Mathews, G. J., & Olive, K. A. 2004a, Phys. Rev. Lett., 92, 
041102 

Ashenfelter, T. P., Mathews, G. J., & Olive, K. A. 2004b, Astrophys. J., 615, 82 
Audi, G, H., W. A., & C, T. 2003, Nucl.Phys. A, 729 

Barrow, J. 2002, The constants of nature: from alpha to omega (London, 

England: Jonathan Cape) 
Bean, S. & Savage, M. 2003, Nucl. Phys. A, 717, 91 
Bergstrom, L., Iguri, S., & Rubinstein, H. 1999, Phys. Rev., D60, 045005 
Brans, C. & Dicke, R. H. 1961, Phys. Rev., 124, 925 

Brink, D. 1966, Proc. Int. School "Enrico Fermi" 36, Varenna 1965, Academic 

Press, New- York, 247 
Bromm, V., Yoshida, N., Hernquist, L., & McKee, C. F. 2009, Nature, 459, 49 
Calmet, X. & Fritzsch, H. 2002, Phys. Lett., B540, 173 
Campbell, B. A. & Olive, K. A. 1995, Phys. Lett. B, 345, 429 
Chand, H., Srianand, R., Petitjean, P., & Aracil, B. 2004, A&A, 417, 853 
Christlieb, N., Gustafsson, B., Korn, A. J., et al. 2004, ApJ, 603, 708 
Coc, A., Nunes, N. J., Olive, K. A., Uzan, J.-R, & Vangioni, E. 2007, 

Phys. Rev. D, 76, 023511 
Cook, C. W., Fowler, W. A., Lauritsen, C. C, & Lauritsen, T. 1957, Phys. Rev., 

107, 508 

Csoto, A., Oberhummer, H., & Schlattl, H. 2001, Nucl. Phys., A688, 560 

Damour, T. & Donoghue, J. F. 2008, Phys. Rev. D, 78, 014014 

Damour, T. & Dyson, F. 1996, Nucl. Phys. B, 480, 37 

Damour, T. & Esposito-Farese, G. 1992, Class. Quant. Grav., 9, 2093 

Damour, T. & Nordtvedt, K. 1993, Phys. Rev. Lett., 70, 2217 

Damour, T., Piazza, F, & Veneziano, G. 2002a, Phys. Rev. Lett., 89, 081601 

Damour, T., Piazza, F, & Veneziano, G. 2002b, Phys. Rev. D, 66, 046007 

Damour, T. & Polyakov, A. M. 1994, Nucl. Phys. B, 423, 532 

Dent, T. & Fairbairn, M. 2003, Nucl. Phys. B, 653, 256 

Dent, T., Stern, S., & Wetterich, C. 2007, Phys. Rev., D76, 063513 

Dicke, R. H. 1959, Nature, 183, 170 

Dine, M., Nir, Y, Raz, G, & Volansky, T. 2003, Phys. Rev., D67, 015009 
Dirac, P. A. M. 1937, Nature, 139, 323 

Dmitriev, V. F. & Flambaum, V. V. 2003, Phys. Rev. D, 67, 063513 
Dmitriev, V. F, Flambaum, V. V., & Webb, J. K. 2004, Phys. Rev. D, 69, 063506 
Duff, M. J. 2002, arXiv:hep-th/0208093 
Duff, M. J., Okun, L. B., & Veneziano, G. 2002, JHEP, 3, 23 
Dunbar, D. N., Pixley, R. E., Wenzel, W. A., & Whaling, W. 1953, Phys. Rev., 
92,649 

Dyson, F. J. 1972, Aspects of Quantum Theory, edited by A. Salam and E.P. 

Wigner (Cambridge U. Press), 213 
Ekstrom, S., Meynet, G, Chiappini, C, Hirschi, R., & Maeder, A. 2008, A&A, 

489, 685 

Epelbaum, E., Meissner, U. G, & Gloeckle, W. 2003, Nucl. Phys. A, 714, 535 
Fairbairn, M. 1999, astro-ph/99 10328 
Fierz, M. 1956, Helv. Phys. Acta, 29, 128 

Flambaum, V. V. & Shuryak, E. V. 2002, Phys. Rev., D65, 103503 
Flambaum, V. V. & Shuryak, E. V. 2003, Phys. Rev. D, 67, 083507 
Frebel, A., Collet, R., Eriksson, K., Christlieb, N., & Aoki, W. 2008, ApJ, 684, 
588 

Fujii, Y. & Iwamoto, A. 2003, Phys. Rev. Lett., 91, 261101 
Fujii, Y. et al. 2000, Nucl. Phys., B573, 377 

Fynbo, H. O. U., Diget, C. A., Bergmann, U. C, et al. 2005, Nature, 433, 136 

Gamow, G. 1967, Phys. Rev. Lett., 19, 759 

Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 

Hirschi, R. 2007, A&A, 461, 571 

Hoyle, F. 1954, ApJS, 1, 121 

Ichikawa, K. & Kawasaki, M. 2002, Phys. Rev., D65, 123511 
Ichikawa, K. & Kawasaki, M. 2004, Phys. Rev, D69, 123506 
Iliadis, C. 2007, Nuclear Physics of Stars (Wiley-VCH Verlag, Wenheim, 
Germany) 

Jordan, P. 1937, Naturwissenschaften, 25, 513 
Jordan, P. 1949, Nature, 164, 637 

Kaluza, T. 1921, Preuss. Akad. Wiss. Berlin (Math. Phys.) K, 1, 966 
Kamimura, M. 1981, Nucl. Phys. A, 351, 456 
Khatri, R. & Wandelt, B. D. 2007, Phys. Rev. Lett., 98, 111301 
Klein, O. 1926, Zeitschrift fur Physik, 37, 895 

Kolb, E. W., Perry, M. J., & Walker, T. P. 1986, Phys. Rev. D, 33, 869 

Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 

Korennov, S. & Descouvemont, P. 2004, Nucl. Phys. A, 740, 249 

Kunz, R., Fey, M., Jaeger, M., et al. 2002, ApJ, 567, 643 

Landau, S. J., Mosquera, M. E., & Vucetich, H. 2006, Astrophys. J., 637, 38 

Langacker, P., Segre, G, & Strassler, M. J. 2002, Phys. Lett. B, 528, 121 



Langanke, K., Wiescher, M., & Thielemann, F.-K. 1986, Zeitschrift fur Physik 
A, 324, 147 

Lee, S., Olive, K. A., & Pospelov, M. 2004, Phys. Rev., D70, 083503 
Livio, M., Hollowell, D., Truran, J. W., & Weiss, A. 1989, Nature, 340, 281 
Marigo, P., Girardi, L., Chiosi, C, & Wood, P. R. 2001, A&A, 371, 152 
Martins, C. J. A. P. et al. 2004, Phys. Lett., B585, 29 
Muller, C. M., Schafer, G, & Wetterich, C. 2004, Phys. Rev., D70, 083504 
Murphy, M. T., Webb, J. K., & Flambaum, V. V. 2003, MNRAS, 345, 609 
Murphy, M. T, Webb, J. K., & Flambaum, V. V. 2007, Phys. Rev. Lett., 99, 
239001 

Nakashima, M., Nagata, R., & Yokoyama, J. 2008, Prog. Theor. Phys., 120, 1207 
Navratil, P., Quaglioni, S., Stetcu, I., & Barrett, B. R. 2009, Journal of Phys. G, 
36,083101 

Nollett, K. M. & Lopez, R. E. 2002, Phys. Rev., D66, 063507 
Nomoto, K, Thielemann, F.-K., & Miyaji, S. 1985, A&A, 149, 239 
Oberhummer, H., Csoto, A., Fairbairn, M., Schlattl, H., & Sharma, M. M. 2003, 

Nucl. Phys., A719, 283 
Oberhummer, H., Csoto, A., & Schlattl, H. 2000, Sci, 289, 88 
Oberhummer, H., Csoto, A., & Schlattl, H. 2001, Nucl. Phys., A689, 269 
Olive, K. A., Pospelov, M., Qian, Y, et al. 2002, Phys. Rev. D, 66, 045022 
Olive, K. A. et al. 2004, Phys. Rev., D69, 027701 

Petrov, Y. V., Nazarov, A. I., Onegin, M. S., Petrov, V. Y, & Sakhnovsky, E. G. 

2006, Phys. Rev. C, 74, 064610 
Pochet, T, Pearson, J., Beaudet, G, & Reeves, H. 1991, Astron. Astrophys., 243, 

1 

Quast, R., Reimers, D., & Levshakov, S. A. 2004, A&A, 415, L7 
Riazuelo, A. & Uzan, J.-P. 2002, Phys. Rev. D, 66, 023525 
Rosenband, T., Hume, D. B., Schmidt, P. O., et al. 2008, Sci, 319, 1808 
Rozental', I. L. 1988, Big bang, big bounce : how particles and fields drive cos- 
mic evoiution (Berlin ; New York : Springer- Verlag, cl988.) 
Schlattl, H, Heger, A., Oberhummer, H, Rauscher, T, & Csoto, A. 2004, 
Ap&SS, 291,27 

Scoccola, C. G, Mosquera, M. E., Landau, S. J., & Vucetich, H. 2008, Phys. 

Rev, D78, 083527 
Shlyakhter, A. f. 1976, Nature, 264, 340 

Srianand, R., Chand, H, Petitjean, P., & Aracil, B. 2004, Phys. Rev. Lett., 92, 
121302 

Srianand, R., Chand, H, Petitjean, P., & Aracil, B. 2007, Physical Review 

Letters, 99, 239002 
Stefanescu, P. 2007, New Astron., f 2, 635 

Suzuki, Y, Matsumura, H, Orabi, M., et al. 2008, Phys. Lett. B, 659, 160 

Taylor, T. R. & Veneziano, G. 1988, Phys. Lett. B, 213, 450 

Thompson, D. R., Lemere, M., & Tang, Y. C. 1977, Nucl. Phys. A, 286, 53 

Tilley, D., Kelleya, J., Godwina, J., et al. 2004, Nucl.Phys. A, 745 

Tur, C, Heger, A., & Austin, S. M. 2007, ApJ, 671, 821 

Uzan, J.-P. 1999, Phys. Rev. D, 59, 123510 

Uzan, J.-P. 2003, Rev. of Mod. Phys., 75, 403 

Uzan, J.-P. 2005, AIP Conf. Proc, 736, 3 

Uzan, J.-P. 2009a, Gen. Ref. Grav., in press, arXiv:astro-ph/0908.2243 

Uzan, J.-P. 2009b, Space Sci. Rev., in press, arXiv:gr-qc/0907.308i 

Uzan, J.-P. & Leclercq, B. 2008, The natural laws of the universe: Understanding 

fundamental constants (Berlin, Germany: Springer) 
Webb, J. K, Murphy, M. T, Flambaum, V. V, et al. 2001, Phys. Rev. Lett., 87, 

091301 

Weinberg, S. 1983, in Phil. Trans. R. Soc. London, Series A, Vol. 310, 249-252 
Wetterich, C. 1988, Nucl. Phys., B302, 645 
Wetterich, C. 2003, Phys. Lett. B, 561, 10 

Wildermuth, K. & Tang, Y. C. 1977, A unified theory of the nucleus (Vieweg, 
Braunschweig) 

Will, C. 1993, Theory and experiment in gravitational physics (Cambridge 

University Press, England) 
Witten, E. 1984, Phys. Lett. B, 149, 351 
Woosley, S. E. & Weaver, T. A. 1986, ARA&A, 24, 205 
Wu, Y.-S. & Wang, Z. 1986, Phys. Rev. Lett., 57, 1978 
Yoo, J. J. & Scherrer, R. J. 2003, Phys. Rev. D, 67, 043517 



